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Abstract. Algorithm, based on explicit representations for the analytic contin- 
uation of Faddeev components of the three-body T- matrix in unphysical energy 
sheets, is employed to study mechanism of disappearance and formation of the 
Efimov levels of the helium 4 He trimer. 

1 Introduction 

Explicit representations for the Faddeev components of the three-body T- 
matrix continued analytically into unphysical sheets of the energy Riemann 
surface have been formulated and proved recently in Ref. Jlj] . According to the 
representations, the T-matrix in unphysical sheets is explicitly expressed in 
terms of its components only taken in the physical sheet. Analogous explicit 
representations were also found for the analytic continuation of the three-body 
scattering matrices. These representations disclose the structure of kernels of 
the T- and S-matrices after continuation and give new capacities for analytical 
and numerical studies of the three-body resonances. In particular the represen- 
tations imply that the resonance poles of the S-matrix as well as T-matrix in 
an unphysical sheet correspond merely to the zeros of the suitably truncated 
three-body scattering matrix taken in the physical sheet. Therefore, one can 
search for resonances in a certain unphysical sheet staying always, nevertheless, 
in the physical sheet and only calculating the position of zeros of the appropri- 
ate truncation of the total three-body scattering matrix. This statement holds 
true not only for the case of the conventional smooth quickly decreasing inter- 
actions but also for the case of the singular interactions described by different 
variants of the Boundary Condition Model, in particular for the inter-particle 
interactions of a hard-core nature like in most molecular systems. 

'Contribution to Proceeding of the 16th European Confe rence on Fcw-RoH y Problems in 
Physics, Autrans (France), 1-6 June 1998. LANL E-print [pnysics/9810004 
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As a concrete application of the method, we present here the results of 
our numerical study of the simplest truncation of the scattering matrix in the 
4 He three- atomic system, namely of the (2 + 1— >2 + l) S-matrix component 
corresponding to the scattering of a 4 He atom off a 4 He dimer. The point is 
that there is already a series of works [§J-@ (also see Refs. §-§) showing that 
the excited state of the 4 He trimer is initiated by the Efimov effect || . In these 
works, various versions of the He- He potential were employed. However, the 
basic result of Refs. on the excited state of the helium trimer is the 

same: this state disappears after the interatomic potential is multiplied by the 
increasing factor A when it approaches the value about 1.2. It is just such 
a nonstandard behavior of the excited-state energy as the coupling between 
helium atoms becomes more and more strengthening, points to the Efimov 
nature of the trimer excited state. The present work is aimed at elucidating 
the fate of the trimer excited state upon its disappearance in the physical 
sheet when A > 1 and at studying the mechanism of arising of new excited 
states when A < 1. As the interatomic He -He potential, we use the potential 
HFD-B H . We have established that for such He - He - interactions the trimer 
excited-state energy merges with the two-body threshold at A « 1.18 and 
with further decreasing A it transforms into a virtual level of the first order (a 
simple real pole of the analytic continuation of the (2 + 1 — > 2 + 1) scattering 
matrix component) lying in the unphysical energy sheet adjoining the physical 
sheet along the spectral interval between ed and the three-body threshold. We 
trace the position of this level for A increasing up to 1.5. Besides, we have 
found that the excited (Efimov) levels for A < 1 also originate from virtual 
levels of the first order that are formed in pairs. Before a pair of virtual levels 
appears, there occurs a fusion of a pair of conjugate resonances of the first order 
(simple complex poles of the analytic continuation of the scattering matrix in 
the unphysical sheet) resulting in the virtual level of the second order. 

2 Representations for three-body T and S matrices 
in unphysical energy sheets 

The method used for calculation of resonances in the present work, is based on 
the explicit representations |l[] for analytic continuation of the T- and scattering 
matrices in unphysical sheets which hold true at least for a part of the three- 
body Riemann surface. To describe this part we introduce the auxiliary vector- 
function f(z) = (f (^),fi,i(^),...,fi, ni (^), h.i(z), ...,f 2 ,„ 2 (z), fs,i(z),-,h,n 3 (z)) 
with fo(z) = lnz and f a ,j(z) = {z — A^j) 1 / 2 . Here, by z we understand the 
total three-body energy in the c. m. system and by A aj -, the respective binding 
energies of the two-body subsystems a, a = 1, 2, 3, j = 1, 2, n a , n a < oo. The 
sheets n; of the Riemann surface of the vector-function f (z) are numerated by 
the multi-index I = (Z , h,i, — , h,ni, h,i, —, h,n 2 ,h,i, —, h,n 3 ), where l a j = if 
the sheet n; corresponds to the main (arithmetic) branch of the square root 
(z — Aq,) 1 / 2 . Otherwise, l a j = 1 is assumed. Value of Iq coincides with the 
number of the branch of the function lnz, lnz = ln|z| + i2ttIq + i</> where 
<j) — argz. For the physical sheet identified by Iq — l a j = 0, a — 1,2,3, j = 
1, 2, n a , we use the notation Ho- 

Surely, the structure of the total three-body Riemann surface is essentially 
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more complicated than that of the auxiliary function f. For instance, the sheets 
II; with Iq = ±1 have additional branching points corresponding to resonances 
of the two-body subsystems. The part of the total three-body Riemann surface 
where the representations of Ref. jy are valid, consists of the sheets II; of the 
Riemann surface of the function f identified by Iq = (such unphysical sheets 
are called two-body sheets) and two three-body sheets identified by = ±1 and 
l a ,j = 1, a = 1,2, 3, j = 1,2, ...,n a . 

In what follows by k a ,p a {k ai p a £ R 3 or k a ,p a S C 3 ) we under- 
stand the standard reduced relative momenta of the three-body system while 
P = {k a ,p a } (P £ R 6 or P € C 6 ) stands for the total relative momentum. 

The representations jl| for the analytic continuation of the matrix M(z) = 
{M a p(z)} , a, j3 = 1,2,3, of the Faddeev components Map(z) (see of the 
three-body T-operator, into the sheet IT; read as followsQ 



M(z) 



= M(z) - B\z)A{z)L^ 1 {z)LB{z). (1) 



l 

n 



Here, the factor A(z) is the diagonal matrix, A(z) — diag{Ao(z), Ai t i{z), . . . , 
Ai.n 3 (z)} with Aq(z) — —iriz 2 and A a j = —iri^/ z — \ a .j, j = 1,2, ...,n a . 
Notations L and L stand for the diagonal number matrices combined of the in- 
dices of the sheet IT;: L = diag{Z , Ji i, h,n 3 } and L = diag{ |Z 1 , l3, n3 }- 
By S;(z) we understand a truncation of the three-body scattering matrix S(z) 

*~ 3 n a 

defined in Q = L 2 (S 5 ) © © L 2 (S 2 ) by the equation 

Q = l j = l 

S l {z)=f+L[S(z)-T\L 
where / is the identity operator in Q . Also, we use the notations 

B(z) = ( j^ XM + v] ) and Bt(z) = (m(z)QUI [v + MT]*jt 

Here, v = diag{vi, v 2 , v^} with v a , the pair potentials, a = 1, 2, 3. At the same 
time, fi = (1, 1, 1), T = {Y a p} with T a p = 1 — 8 a p-, a,f3= 1,2,3, and 
= diaglM/i, S&2, ^3} where \& Q , a = 1,2,3, are operators acting on / = 

(/1,/2,-,/nJ € © L 2 (R 3 ) as (* a /)(P) = J2 i>a,j(k a )fj(p a ) where, in turn, 
j=i j=i 

ip a j is the bound-state wave function of the pair subsystem a corresponding 
to the binding energy \ a ,j- By Vl/* we denote operator adjoint to \I>. Notation 
Jo(z) is used for the operator restricting a function on the energy-shell \P\ 2 = z. 
The diagonal matrix-valued function Ji(z) — diag{Ji.i(z), J3,n 3 (^)} consists 
of the operators J Q ,j(z) of restriction on the energy surfaces \p a \ 2 = z — X a ,j- 
The operators fit, 3q(z) and jj (z) represent the "transposed" matrices fi, Jq(z) 
and Ji(z), respectively. Operators Jq(z) and J{(z) act in the expression for E$t 
(as if) to the left. 

1 One assumes that all the pair interactions fall off in the coordinate space not slower than 
exponentially and, thus, their Fourier transforms v a (k), k S C 3 , are holomorphic functions 
of the relative momenta k in a stripe | Im k\ < b for some b > 0. 
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With some stipulations (see jjj) the representations for the scattering matrix 
read 



S(z) 



= £ 



(l){T+S- 1 (z)[S(z)~T}e(l)}£(l). 



(2) 



Here, £ = diag{£Q)£i,i) ...,£3^3} where £q is the identity operator in 1/2(5' ) if 
lo = and inversion, (£ f)(P) = f(—P), if lo = ±1- Analogously, £ a j is the 
identity operator in L2{S 2 ) for l a j = and inversion for l a j = 1. Notation 
e(l) is used for the diagonal number matrix e(l) = diagjeo, ex,i, ...,e3.„ 3 } with 
nontrivial elements e QJ = 1 if i a j = and e a j — — 1 if Z Q j = 1; for all the 
cases eo = 1. 

It follows from the representations ([!]) and (||) that the resonances (the 



nontrivial poles of M(z) 
are those points z = 



and S(z) 



) situated in the unphysical sheet II; 



rcs in the physical sheet where the matrix S/(z) has zero 
as eigenvalue. Therefore, calculation of resonances in the unphysical sheet Hi is 
reduced to a search for zeros of the respective truncation Sj(z) of the scattering 
matrix S(z) in the physical sheet. 



3 Method for search of resonances in a three body system 
on the basis of the Faddeev differential equations 

In this work we discuss the example of the three-atomic 4 He system at the 
total angular momentum L = 0. We consider the case where the interatomic 
interactions include a hard core component and, outside the hard core domain, 
are described by conventional smooth potentials. In this case, the angular par- 
tial analysis reduces the initial Faddeev equation for three identical bosons to 
a system of coupled two-dimensional integro-differential equations (see Rcf. ^ 
and references therein) 



d 2 d 2 ( 1 1 

9a; 2 dy 2 \x 2 y 2 



Fi(x,y) (3) 



-V(x)Wi(x,y), x>c 
0, x < c . 

Here, x, y stand for the standard Jacobi variables and c, for the core range. 
At L = the partial angular momentum / corresponds both to the dimcr 
subsystem and a complementary atom. For the iS-state three-boson system I is 
even, I = 0, 2, 4, . . . . In our work, the energy z can get both real and complex 
values. The He-He potential V(x) acting outside the core domain is assumed 
to be central. The partial wave function <Pi(x,y) is related to the Faddeev 
components F t (x,y) by &t(x,y) = F t (x,y) + /+ drj h W (x, y, v) Fv v') 

where x' = (± x 2 + § y 2 - & xyrj) 1 / 2 ,y' = (^x 2 + \ y 2 + & X yn) 1/2 and 
1 < rj < 1. The explicit form of the function hw can be found in Refs. [[Io| Ell. 
The functions Fi(x,y) satisfy the boundary conditions 

Fi(x,y)\ x=0 =F l (x,y)\ y=0 =0, iP l (c,y) = (4) 
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Note that the last of these conditions is a specific condition corresponding to 
the hard-core model (see Ref. [||). 

Here we only deal with a finite number of equations (^), assuming that 
I < (max where Z max is a certain fixed even number. The condition < / < i max 
is equivalent to the supposition that the potential V(x) only acts in the two- 
body states with I = 0,2, ... , l ma , x . The spectrum of the Schrodingcr operator 
for a system of three identical bosons with such a potential is denoted by o^b- 
We assume that the potential V(x) falls off exponentially and, thus, |V^(a;)| < 
C exp(— fix) with some positive C and fi. For the sake of simplicity we even 
assume sometimes that V(x) is finite, i. e., V(x) = for x > r*o, ro > 0. Looking 
ahead, we note that, in fact, in our numerical computations of the 4 He3 system 
at complex energies we make a "cutoff" of the interatomic He - He - potential 
at a sufficiently large tq. The asymptotic conditions as p — > oo and/or y — ► oo 
for the partial Faddeev components of the (2 + 1 — > 2 + 1 ; 1 + 1 + 1) scattering 
wave functions for z = E + iO, E > 0, read (see, e. g., Ref. JlCfl) 

Fi{x,y;z) = 5i ip d (x) {sin(Vz - e d y) + cxp(Vz - e d y) [&o(z) +o(l)]} 



We assume that the 4 He dimer has an only bound state with an energy e d , 
e d < 0, and wave function tp d (x), ipd(x) — for < x < c. The notations p, 
p = \J x 2 + y 2 , and 8, 8 = a,rct&n(y/x), are used for the hyperradius and hy- 
perangle. The coefficient a.o(z), z = E + iO, for E > e d is the elastic scattering 
amplitude. The functions Ai(E + iO, 8) provide us, at E > 0, the correspond- 
ing partial Faddeev breakup amplitudes. For real z = E + iO, E > e d , the 
(2 + 1— >2 + 1) component of the s-wave partial scattering matrix for a system 
of three helium atoms is given by the expression 



Our goal is to study the analytic continuation of the function So(^) into the 
physical sheet. As it follows from the results of Refs. [|), the function So (z) is 
just that truncation of the total scattering matrix whose roots in the physical 
sheet of the energy z plane correspond to the location of resonances situated 
in the unphysical sheet adjoining the physical one along the spectral interval 



There are the following three important domains in the the physical sheet. 

1°. The domain 11^ where the Faddeev components Fi(x, y; z) (and, hence, 
the wave functions &i(x, y; z)) can be analytically continued in z so that the 
differences $i(x,y;z) = Fi(x,y;z) - Si Q ip d (x) sin(^z - e d y) at z £ iT (<fr) \<r 3B 
are square integrable. This domain is described by the inequality 



+ 



exp(i A /i / o) 



[A l (z,e) + o(i)] . 



(5) 



S (z) = l + 2ia (z). 



(cd,0). 



Im \fz — e d < min 




2°. The domain 11^ where both the elastic scattering amplitude a (z) and 
the Faddeev breakup amplitudes Ai(z,9) can be analytically continued in z, 
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z ^ 03B, and where the continued functions i*z(x, y; z) still obey the asymptotic 
formulas This domain is described by the inequalities 

Im \[z + Im y/ z - e d < — p . 

3°. And finally, we distinguish the domain IT^ S ', most interesting for us, 
where the analytic continuation in z, z ^ c^b, can be only done for the am- 
plitude ao(z) (and consequently, for the scattering matrix So(z)); the analytic 
continuabilty of the amplitudes Ai(z, 6) in the whole domain II^ S ' is not re- 
quired. The set II^ S ' is a geometric locus of points obeying the inequality 



Im ^Jz - e d < min < —=■ yj\e d \, — — 




Since the spherical wave exp{\yfz p)j ' ^fp in Eq. (|J) is a function rapidly de- 
creasing in all the directions, the use of the asymptotic condition (|^) is justified 
even if z 6 11^ \H { - AS > . Outside of the domain the numerical construction 
of So (z) by solving the Faddeev differential equations is, in general, impossible. 



4 Numerical results 



In the present work we search for the resonances of the 4 He trimer including 
the virtual levels as roots of So(z) and for the bound-state energies as positions 
of poles of So(z). All the results presented below are obtained for the case 
'max = 0. In all our calculations, h 2 /m = 12.12 K A 2 . As the interatomic He- 
He - interaction we employed the HFD-B potential constructed by R. A. Aziz 
and co-workers ||. 

The value of the core range c is chosen to be so small that its further 
decrease does not appreciably influence the dimer binding energy ed and the 
trimer ground-state energy e[ 0) . Unlike the paper W], where c was taken to 
be equal 0.7 A, now we take c = 1.3 A. We have found that such a value of c 
provides at least six reliable figures of e d and three figures of . 

Since the statements of Sect. || are valid, generally speaking, only for the 
potentials decreasing not slower than exponentially, we cut off the potential 
HFD-B setting V(x) = for x > r$. We have established that this cutoff for 

r > 95A provides the same values of (ed — — 1.68541mK), = 
— 0.096 K) and scattering phases which were obtained in earlier calculations 
[|| performed with the potential HFD-B. Also, we have found that the trimer 
excited-state energy E^ = — 2.46mK. Comparison of these results with results 
of other researchers can be found in Ref . Q . In all the calculations of the present 
work we take ro = 100 A. Note that if the formulas from Sect. || describing the 
holomorphy domains II^ A ' and TT^- 1 are used for finite potentials, one 

should set in them /! = +00. 

A detailed description of the numerical method we use is presented in 
Ref. [Q . When solving the boundary- value problem we carry out its finite- 
difference approximation in polar coordinates p and 9. In this work, we used 
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Figure 1. Root locus curves of the real and imaginary parts of the scattering matrix 
So (2). The solid lines correspond to Re So (2) = 0, while the tiny dashed lines, to 
ImSo(z) = 0. The numbers 1, 2 and 3 denote the boundaries of the domains 
II^ S ' and II^ A \ respectively. Complex roots of the function So (2) are represented by 
the intersection points of the curves Re So(z) — and Im So (2) = and are located at 
(-2.34 + i0.96) mK, (-0.59 + i2.67) mK, (2.51 + i4.34) mK and (6.92 + i6.10) mK. 



the grids of dimension N§ — N p — 600 — 1000. In essential, we chose the values 
of the cutoff hyperradius p max = Pn p from the scaling considerations (sec [|j ) . 
We solve the resultant block-three-diagonal algebraic system on the basis of 
the matrix sweep method. This allows us to dispense with writing the system 
matrix on the hard drive and to carry out all the operations related to its inver- 
sion immediately in RAM. Besides, the matrix sweep method reduces almost 
by one order the computer time required for computations on the grids of the 
same dimensions as in 0]. 

Because of the symmetry relationship So{z) — Sq(~z) we performed all the 
calculations for So(z) only at Imz > 0. First, we calculated the root lines 
of the functions ReSo(z) and ImSo(.2). For the case of the grid parameters 
N$ = N p = 600 and p max = 600 A these lines are depicted in Fig. 1. Both res- 
onances (roots of S (z)) and bound-state energies (poles of S (z)) of the 4 He 
trimer are associated with the intersection points of the curves Re So(z) = and 
ImSo(z) = 0. In Fig. 1, along with the root lines we also plot the boundaries of 
the domains II^ S \ n {A ~> and 77^. One can observe that the "good" domain 
Tj(-S) i nc l u des none of the points of intersection of the root lines Re So (z) = 
and ImSo(z) = 0. The caption for Fig. 1 points out positions of four "res- 
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Table 4.1. Dependence of the dimer binding energy ed and the differences ed — E\ , 
ed — E\ , €d — E t and — E\ (all in mK) between this energy and the trimer 
exited-state energies E^ , E^ and the virtual-state energies E^* , E^** on the 
factor A. 



A 


ed 






e d - EF>~ 




0.995 


-1.160 


0.710 








0.990 


-0.732 


0.622 








0.9875 


-0.555 


0.573 


0.473 


0.222 




0.985 


-0.402 


0.518 


0.4925 


0.097 




0.980 


-0.170 


0.39616 


0.39562 


0.009435 




0.975 


-0.036 


0.2593674545 


0.2593674502 




0.00156 



onances", the roots of So(z), found immediately beyond the boundary of the 
domain ]j( s \ It is remarkable that the "true" (i. e., getting inside II^ ) virtual 
levels and then the energies of the excited (Efimov) states appear just due to 
these (quasi)resonances when the potential V{x) is weakened. 

Following |2|]-j§], instead of the initial potential V(x) = Vhfd-b(£), we, 
further, consider the potentials V{x) = A ■ Vhfd-b(^)- To establish the mech- 
anism of formation of new excited states in the 4 He trimer, we have first cal- 
culate the scattering matrix So(z) for A < 1. In Table 4.1 for some values of 
A from the interval between 0.995 and 0.975, we present the positions of roots 
and poles of So(z), we have obtained at real z < £d(A). We have found that 
for a value of A slightly smaller than 0.9885, the (quasi)resonance closest to 
the real axis (see Fig. 1) gets on it and transforms into a virtual level (the 
root of So(z)) of the second order corresponding to the energy value where the 
graph of So(z), z G K, z < ed, is tangent to the axis z. This virtual level is 
preceded by the (quasi)resonances z = (—1.04 + i 0.11) mK for A = 0.989 and 
z = (-0.99 + i0.04) mK for A = 0.9885. With a subsequent decrease of A the 
virtual level of the second order splits into a pair of the first order virtual levels 
E^* and E^** , E^* < E^** which move in opposite directions. A charac- 
teristic behavior of the scattering matrix So (z) when the resonances transform 

into virtual levels is shown in Fig. 2. The virtual level E t moves towards 
the threshold ed and "collides" with it at A < 0.98. For A = 0.975 the function 
Sq(z) has no longer the root corresponding to E^** . Instead of the root, it 
acquires a new pole corresponding to the second excited state of the trimer 

(2) 

with the energy E t . We expect that the subsequent Efimov levels originate 

from the virtual levels just according to the same scheme as the level E). ' does. 

The other purpose of the present investigation is to determine the mecha- 
nism of disappearance of the excited state of the helium trimer when the two- 
body interactions become stronger owing to the increasing A > 1. It turned out 
that this disappearance proceeds just according to the scheme of the formation 
of new excited states; only the order of occurring events is inverse. The results 
of our computations of the energy E^ when A changes from 1.05 to 1.17 are 



8.3 




Figure 2. Graphs of the function So(z) at real z < ed for three values of A < 1 . The 
notations used: E* = E? ] * /\e d \, E" = E 4 (2) *7|e d |. 



Table 4.2. Dependence of the dimer energy e^, the difference ed — E^ between ed 
and the trimer exited-state energy E\ and the difference ed — E] between e<j and 
the trimer virtual-state energy E^ (all in mK) on the factor A. 



A 




£d - E t 


A 






1.05 


-12.244 


0.873 


1.18 


-82.927 


0.001 


1.10 


-32.222 


0.450 


1.20 


-99.068 


0.057 


1.15 


-61.280 


0.078 


1.25 


-145.240 


0.588 


1.16 


-68.150 


0.028 


1.35 


-261.393 


3.602 


1.17 


-75.367 


0.006 


1.50 


-490.479 


12.276 



given in Table 4.2. In the interval between A = 1.17 and A = 1.18 there occurs 
a "jump" of the level on the unphysical sheet and it transforms from the 

pole of the function So(z) into its root, E\ ? , corresponding to the trimer vir- 
tual level. The results of calculation of this virtual level where A changes from 
1.18 to 1.5 are also presented in Table 4.2. 

More details of our techniques and material presented will be given in an 
extended article Q. 
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